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Abstract 

By expressing thermodynamic functions in terms of tlie edge and density of Lee- Yang ze- 
roes, we relate the scaling behaviour of the specific heat to that of the zero field magnetic 
susceptibility in the thermodynamic limit of the Xy-model in two dimensions. Assuming that 
finite-size scaling holds, we show that the conventional Kosterlitz-Thouless scaling predictions 
for these thermodynamic functions are not mutually compatable unless they are modified by 
multiplicative logarithmic corrections. We identify these logarithmic corrections analytically in 
the case of the specific heat and numerically in the case of the susceptibility. The techniques 
presented here are general and can be used to check the compatibility of scaling behaviour of 
odd and even thermodynamic functions in other models too. 
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The Kosterlitz— Thouless Scenario The partition function for the 0{n) non-Unear a-model 
defined on a regular d-dimensional lattice A with periodic boundary conditions in the presence of 
an external magnetic field can be written as 



{CTx} 

where n is a unit vector defining the direction of the external magnetic field and /i is a scalar 
parameter representing its strength. The summation is over all configurations open to the system. 
The factor J\f ensures that Zl(0,0) = 1 and 

d 

xeAfi=i xeA 

Here L represents the linear extent of the system, /? = 1/kT and ffx is an n-component unit length 
spin at site x £ A. The case d = 2, n = 2 is the two dimensional Xy-model or plane rotator 
model and is unusual in that it exhibits an exponentional singularity. Using an approximate 
renormalization group approach, Kosterlitz and Thouless Q demonstrated the existence of a phase 
transition driven by the condensation of vortices. The model remains critical (thermodynamic 
functions diverge) for all (3 > (3c and the critical exponents are dependent on temperature. In 
terms of the reduced temperature 

the scaling behaviour of the correlation length, susceptibility and the specific heat is given in as 

ioo{t) ~ e^'-\ (2) 
Xoo{t) ~ (3) 
Coo(i) ~ C + constant, (4) 

where for t— >0~'',z^ = l/2,ry = l/4 and a = —d = —2. The purpose of this work is to argue that the 
latter two scaling formulae are incompatible as they stand. To this end, a method is presented by 
which odd and even thermodynamic functions (the susceptibility and specific heat) can be related 
and expressed in terms of partition function zeroes. Using certain reasonable assumptions regarding 
finite-size scaling, it is shown that there have to exist multiplicative logarithmic corrections ^ to 
d) and (|). This method can be applied to any model. 



Partition Function Zeroes Lee and Yang showed the connection between critical behaviour 

and partition function zeroes. For a finite system these are strictly complex (non-real). As L ^ oo 

one expects the zeroes to condense onto smooth curves. Zeroes in the plane of complex external 

magnetic field h are refered to as Lee- Yang zeroes. Lee and Yang further showed that for certain 

systems these zeroes are in fact restricted to the imaginary h axis (the Lee- Yang theorem) |^]. In 

the symmetric phase, t > 0, they lie away from the real /i-axis, pinching it only as t — > 0^ (in the 

^The work of |l[ contains implicitly a prediction for logarithmic corrections but this seems not to have been 
followed up quantitatively. 



thermodynamic limit). For systems obeying the Lee- Yang theorem, such as the XY Model |^], this 
pinching occurs at /i = prohibiting analytic continuation from Re(/i) < to Re(/i) > 0. This 
means that (for the finite-size system or in the thermodynamic limit) the partition function is an 
entire function of h all of whose zeroes are purely imaginary. The zeroes in the complex fugacity 
plane {z = exp (h)) are distributed on the unit circle. Writing these zeroes as 

the partition function is 

ZL{P,h)=pL{P,h)ll{z-e'<>^^'''>), (5) 
j 

where pL is a non- vanishing function of h related to the spectral density (ignored in what follows 
since it contributes only to the regular part of the free energy in the thermodynamic limit). The 
free energy is then 

/L(/3,/i)=L-'^^ln(z-e*^^(^)). (6) 
j 

Define the density of zeroes to be |5| 

gUP,e) = L-'^Y.^ie - em = • (7) 

j 

Here Gl{(3,0) is the cumulative density of zeroes and monotonically increases from Gl{(3,0) = 
to Gl(/5, vr) = 1/2. The distribution of zeroes is symmetric about the real h- (or z-) axis and so 

gL{P,-0)=gL{P,e). (8) 

In terms of the density of zeroes, (^) becomes 

/^(/3,/t) = i(/i + ln2) + r In {cosh h- cos 9)dGL{(3,e) , (9) 
2 Jo 

where (^ has been used. 

In the high temperature phase there exists a region around 9 = which is free from Lee- Yang 
zeroes Q. We define the corresponding Yang-Lee edge 6'yl(/5) by 

g{p,9) = for - < 9 < 9yM . 

The integral in (P) can therefore be taken to begin at ^yl(/3)- Salmhofer has proved the existence 
of a unique density of zeroes in the thermodynamic limit [^. In this limit, integrating (^ by parts 
and expanding the trigonometric functions gives for the singular part of the free energy ^, ^ 



/oo(/5, h) = -2 ^^Goo(/3, 9)d9 . (10) 

The magnetic susceptibility Xoo is given by the second derivative of the free energy with respect to 
h. Following |, |, |7), this leads to 

G^{f3,9) = Xoo(/3)(^yl(/3))'$ {-Q^y) ' (11) 



^{x) being some function of x such that ^{\ x [< 1) = 0. From ( [10| ) and the fact that G {Oyl, P) = 0, 
one gets the specific heat 



■"'C-C-^dO, (12) 



h=0 

The above formulae are quite general and hold for any model provided it obeys the Lee- Yang 
theorem. 

To proceed further we insert the (model-specific) critical behaviour. Instead of @ and assume 
now the following modified critical behaviour for the singular parts of the zero field susceptibility 
and the specific heat for i > 0. 

Xoo-fa^'^f, C^-et'. (13) 

Similarly, assume that the leading scaling behaviour of the Yang-Lee edge in terms of the reduced 
temperature is 

eYL(t)~eV. (14) 

For t ^ the critical indices are independent of t [0. The scaling behaviour of Coo can be related 
to that of Xoo and Oyl through ( [TlD and (p^). These give 

iitft'i OC ^(i)2-^+2A^2p+r-2-2,. _ 

Therefore, 

A = i(a-2 + r?), (15) 
p = ^(g-r) + l + i^. (16) 

Finite— Size Scaling The finite-size scaling (FSS) hypothesis, first formulated in 1972 by Fisher 
[^, is a relationship between the scaling behaviour of thermodynamic quantities in the infinite 
volume limit and the size dependence of their finite volume counterparts. The general statement of 
FSS, which is expected to hold in all dimensions including the upper critical one g is that if PL(t) 
is the value of some thermodynamic quantity P at reduced temperature t measured on a system of 
linear extent L, then 

^l(O) _ / a(0) \ 

where ^L{t) is the correlation length of the finite-size system. Luck has shown that, for the XY- 



(17) 



model in two dimensions, ^l(O) is proportional to L [IC]. Fixing the scaling variable, one has 
and from (|^) 

t ~ (lnL)~~ , 

for large enough L. Therefore FSS for the susceptibility and the Yang-Lee edge is 

Xl(0) ~ L^-\\nL)-i, (18) 
^1(0) ~ L^(lnL)-^. (19) 



For the finite size system it is convenient to consider the zeroes in the complex /i-plane. The 
singular part of the free energy corresponds to 



fL{t,h)=L-<'Y,ln{h-ie,{t)) (20) 

j 

where 9i = 9yl- The (reduced) magnetic susceptibility is therefore 

= -L'T.i- (21) 



XL{t) 



dh? 

The lowest lying zeroes are expected to scale in the same way [0] so that 

XL(t)cx-rti(t)-2. (22) 



Prom (|18|) and ([L9|) we have 



and conclude 



L2-''(lnL)~^ ~ L~'^~^^{\nLf-^ , (23) 



a = -d, q=-2{l + iy). (24) 



Thus the scaling behaviour of the singular part of the specific heat indeed exhibits multiplicative 
logarithmic corrections. Accepting the KT predictions u = 1/2, rj = 1/4 for t ^ 0, the leading 
critical exponent for the Yang-Lee edge A and the correction exponents q and p are 

A = -f, g = -3, p = -^-. (25) 

The above analytic considerations have yielded no information on the odd correction exponent 
r. The original renormalisation group analysis of Kosterlitz and Thouless in fact, implicitly 
contained the prediction 

r = -l/16 (26) 



as noted by |12| | . Subsequent analysis have concentrated on the r = form of the scaling behaviour 
(ITsI) and the verification that 7/(/3c) = 1/4. Allton and Hamer |l^ have conjectured that the 
deviation of their determination of rj from 1/4 might be due to logarithmic corrections. 

The study of the full scaling form and the numerical determination of r is the subject of what 
follows. 



Numerical Procedures Numerical methods were used to test the scaling picture of the last sec- 
tion. Specifically, we analysed the FSS behaviour of the Yang-Lee edge so as to test the prediction 

(ID, (H) 

^l(/3e)~L^(lnL)^ (27) 

In particular, we sought to confirm A = 15/8 for the leading behaviour and to check if r 7^ 0. 

An algorithm based on that of Wolff [|l^ was used to simulate the XY-model at zero magnetic 
field {h = 0) on square lattices of sizes L = 32, 64, 128 and 256 . The values of /? at which the 
simulations took place (/3o) were evenly spaced in steps of 0.02 between 1.00 and 1.20. At each 



simulation point, 50, 000 measurements were made from each of two separate starts. The second 
runs were used only to verify equilibration and to improve our understanding of the statistical 
errors. The latter were estimated using bootstrap as described below. The autocorrelation time of 
the susceptibility was also monitored. Histogram reweighting |15] was used to enable extrapolation 
between Po values. 

The determination of Pc is inextricably bound up with the assumed critical scaling behaviour. For 
the present purposes, we adopted a two-stage strategy: we used estimates of f3c from our own 
and other analyses which had assumed no logarithmic corrections then tested the stability of these 
conclusions with respect both to uncertainties in Pc and to the form of the scaling asumption. 

One simple estimate was based on assumed conventional (r = 0) scaling of the susceptibility (P). A 
generalised [^6| Roomany-Wyld approximant |17| was used to give estimates of the corresponding 



renormalisation group beta function for each pair of lattice sizes related by factor of two scaling: 
gRW ^ . 2 r? - 2 + 2r In(ls^)/ In 2 + Hx2l/xl)/ In 2 

A preliminary estimate, using rj = 1/4 and r = 0, leads to to Pc = 1-11(1) which is at the lower 
end of the range of published estimates [jl2|, [|l8|, |9|. These span 1.11 to 1.13. The most recent 
high precision result is 1.1197(5) |]l9[. 



Determination of Lee— Yang Zeroes When the external field is complex {h = hr + ihi), the 
partition function (pi) can be rewritten 



where 



Zl{P, hr + ihi) = ReZi(/3, hr + ihi) + HmZi^P, hr + ihi) , (29) 



ReZiiP, hr + ihi) = jpY. 6^^+'^'-*^ cos (hiM) , 

= Zl{P, hr){cOS (hiM)) fs,hr, (30) 



and 



ImZUP, hr + ihi) = JfY. e^^+^'-^^ sin (hiM) , 

= ZL{P,hr){sm{hiM))p,h,.. (31) 

According to the Lee- Yang theorem the zeroes are on the imaginary /i-axis {hr = 0) where 
{sin (hiM)) j3^hr vanishes and so the Lee- Yang zeroes are simply the zeroes of 

ReZL{P, ihi) oc (cos (/iiM))^,o • (32) 

Thus the Lee- Yang zeroes are easily found. Moreover, one can find them at any P using reweight- 
ing techniques |l^]. For the lowest lying zeroes at P = 1.11 we found 6i{L) = 0.0023353(7), 
0.0006350(2), 0.00017278(5) and 0.000047062(13) for L = 32,64,128 and 256 respectively. Er- 
rors were calculated using the bootstrap method where the data for each Po are resampled 100 



times (with replacement) leading to 100 estimates for 9i, from which the variance and bias can be 
calculated. 

In figure [l| we plot the logarithm of the position of the first Lee- Yang zero against the logarithm 
of the lattice size L, at /? = 1.11. In the absence of any corrections, the slope should give the 
leading power-law exponent A. In fact the slope is -1.8776(2), the deviation from the KT value of 
— 15/8 = —1.875 being attributed to the presence of logarithmic corrections. To identify these, and 
the correction exponent r in (|27|), In (0iL^^/^) is plotted against InlnL in figure |2[ A straight line is 
identified. Its slope is —0.012(1) giving evidence for a non-zero value of r, albeit not in agreement 
with the RG predictions of —1/16 = —0.0625 from [|l|. 

At this point, one must investigate the extent to which these conclusions regarding the existence of 
multiplicative logarithmic corrections depend upon the measurement of (3c- We have attempted to 
do this systematically by studying the above InlnL fits for r as a function of the assumed critical 
beta. Fits with a range of values of r (including zero) are possible but not all with acceptable 
X^/degree of freedom. To be conservative, we take 'acceptable' to mean less than 2. This corre- 
sponds to a minimum confidence level of 14%. In figure ^ we show that acceptable values of 
are possible only for 1.110 ^ /9c ^ 1.120. The r = solution would correspond to /3c ~ 1.105 
and /dof ~ 4. Similary, we find that a fit with r = —1/16 would correspond to j3c ~ 1.138 and 
xVdof « 16. 

In summmary, assumming the Kosteritz-Thouless value r] = 1/4, we find non-zero logarithmic 
corrections to scaling and a corresponding estimate of the critical temperature: 

r = -0.023 ±0.010, /?c = 1.115 ± .005 . (33) 

As a cross check, one may use the above values for r and r] to find the Roomany-Wyld approximant 
(|28|) and estimate (5c from its zeroes. We find (3c = 1.115 it 0.004 in good agreement with the above. 
However, we emphasize the much higher precision available to an analysis based on Lee- Yang zeroes 
rather than on the spin susceptibility. 

Conclusions Theoretical arguments concerning the consistency of the scaling behaviour of odd 
and even thermodynamic functions at a KT phase transition have been presented. The generally 
used scaling formulae have to be modified by multiplicative corrections. These are identified ana- 
lytically for the specific heat and numerically for the susceptibility. This numerical identification 
comes via an analysis of Lee- Yang zeroes, the FSS of which is linked to that of the susceptibility. 

We would like to thank P. Lacock for assistance with multihistogram reweighting techniques. 
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Figure 1: Leading FSS of Lee- Yang Zeroes at (3 = 1.11. 
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Figure 2: Corrections to FSS of Lee- Yang Zeroes at (3 = 1.11. 
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Figure 3: (a) Correction Exponent r Measured over a Range of (3 and (b) corresponding per 
Degree of Freedom (Goodness of Fit). 



